An exactly solvable model for driven dissipative systems 
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We introduce a solvable stochastic model inspired by granular gases for driven dissipative systems. 
We characterize far from equilibrium steady states of such systems through the non-Boltzmann en- 
ergy distribution and compare different measures of effective temperatures. As an example we 
demonstrate that fluctuation-dissipation relations hold, however with an effective temperature dif- 
fering from the effective temperature defined from the average energy. 

PACS numbers: 05.70.Ln, 02.50.Ey, 45.70.-n 



Dissipative many-particle systems are far from thermo- 
dynamic equilibrium, and a general theoretical descrip- 
tion of their statistical mechanics is lacking, in contrast 
to systems in equilibrium, for which there is a well estab- 
lished theory. In this Letter we propose a simple exactly 
solvable model, which provides a context in which certain 
questions concerning driven dissipative systems may be 
resolved unambiguously, such as whether different pro- 
posed "definitions" of temperature give the same value, 
as is the case for thermal equilibrium. 

Granular materials have in recent years been consid- 
ered a paradigm for dissipative open systems . In such 
systems comprised of macroscopic particles, energy is dis- 
sipated via interactions, being transferred from macro- 
scopic degrees of freedom (motion of grains) into micro- 
scopic degrees of freedom, and can not be transformed 
back. Continuous driving is needed in order to maintain 
such a system in a dynamic state. This driving may be 
realized, for instance, by gravity- driven flow down an in- 
cline , by continuous avalanches in a rotating drum Q , 
or in a vibrated container 

An oft-studied system is a homogeneously heated gran- 
ular gas : a collection of hard spheres (or disks in 2D) 
undergoing inelastic collisions and driven by a stochastic 
thermostat. This yields Langevin dynamics including a 
random uncorrelated force F(t), satisfying (F) = and 
(F i (t)F j (t')) = 2~fT B 6{t - t')Sij, and a drag force -jrf. 
This may be interpreted as coupling to a heat bath of 
temperature Tg with a coupling strength 7. Were the 
system non-dissipative, it would reach equilibrium with 
temperature Tb, irrespective of the details of the cou- 
pling, that is, independent of 7. 

A driven dissipative system reaches a non-equilibrium 
steady state, for which a granular temperature Tq may 
be defined as the average kinetic energy of the grains. 
Not only is Tq always smaller than Tb, but its value de- 
pends on the details of the coupling with the bath 0, Q ■ 
This steady state behaves statistically differently than an 
equilibrium state at the same effective temperature. The 
energy distribution deviates from the Boltzmann distri- 
bution, exhibiting overpopulated high energy tails 0. 
Surprisingly, even though driven dissipative systems are 
not in equilibrium, fluctuation-dissipation (FD) relations 



often hold and may serve to define an effective tempera- 
ture, Tfd 0- which has been found in numerical exper- 
iments to coincide with Tq ■ 

Granular gases may be described theoretically by the 
Boltzmann equation of kinetic theory, however solutions 
exist only in terms of approximations valid for small de- 
viations from the equilibrium Boltzmann distribution, 
that is for small inelasticity and low volume fraction pj. 
A dissipative model for which some exact results have 
been found is the one-dimensional Maxwell-model 
in which particles collide inelastically with a uniform col- 
lision rate, independent of velocity or location. Energy 
distributions have also been investigated numerically in 
a two-dimensional version of this model assuming a ran- 
dom impact parameter in every collision . 

In this Letter we present a novel exactly solvable dis- 
sipative model, in which interactions occur randomly 
and redistribute energy randomly between the interact- 
ing particles. Our main results are: 

• All moments of the model may be computed exactly. 

• In the maximally dissipative limit, the generating 
function of the model may be solved exactly. 

• FD relations hold with T FD >T G . 

Our model consists of a collection of N particles hav- 
ing energies Ei, with a constant interaction rate between 
any two particles in the system. In every interaction two 
particles from the system are chosen at random and their 
energies are summed. In the case of conservative dynam- 
ics (analogous to elastic collisions) this total energy is 
repartitioned randomly with a uniform distribution be- 
tween the two interacting particles [12| , while for dissipa- 
tive dynamics (as for inelastic collisions with a constant 
restitution coefficient) only a fraction a of the total en- 
ergy is repartitioned between the particles and the rest is 
dissipated out of the system. Additionally, we couple the 
system to a heat bath so that it reaches a steady state. 

The simplicity of this model results from the fact that 
every particle in it is described only by its energy, as 
opposed, for example, to the 2d degrees of freedom per 
particle in a <i-dimensional frictionless hard sphere gas. 
By eliminating the momentum and spatial variables and 
using solely the energy we turn the vectorial collisions 
between particles into scalar interactions, and preclude 
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spatial correlations. 

For conservative dynamics (a = 1) the system 
reaches an equilibrium state with the exponential Boltz- 
mann distribution for the particles' energies, p(E) = 
T~ 1 exp(—E /T), where the temperature (measured in 
units of energy) equals the average energy in the sys- 
tem T = (E). Dissipative dynamics (a < 1) cause en- 
ergy to decay, and in order to maintain the system in 
a steady state, we keep it in contact with a heat bath, 
which is constructed as an infinitely large system of parti- 
cles obeying the conservative dynamics described above, 
kept in equilibrium at a temperature Tb- The coupling 
of the dissipative system to the bath is through conser- 
vative interactions between a particle chosen at random 
from the system and a particle chosen at random from 
the bath, and is characterized by a coupling strength, /, 
which is defined as the fraction of the particle's inter- 
action that are with the bath out of all its interactions 
(with the bath and with other particles in the dissipative 
system). We will see that the steady state of the dissipa- 
tive system depends on the bath through both Tb and / 

0. 

The stochastic equation of motion for the evolution of 
the energy of particle i during an infinitesimal time step 
dt is hence given by 

{ value : probability : 

Ei(t) 1 - Tdt m 

zaiE^t) + E 3 {t)) (l-f)Tdt ' [1) 
z{Ei{t)+E B ) fTdt 

where: T is the interaction rate per particle per unit 
time (which does not affect the steady state but only the 
rate of approach to it); j 6 {1,...,N} is the index of 
the particle with which particle i may interact, chosen 
randomly at every interaction; z G [0, 1] is the fraction 
of repartitioned energy given to particle i in the inter- 
action, chosen randomly with a uniform distribution at 
every interaction; Eb is the energy of the particle from 
the bath with which particle i may interact, which at ev- 
ery interaction is chosen randomly from the equilibrium 
distribution in the bath, p B {E B ) = T B 1 exp{-E B /T B ). 

We first demonstrate how all moments of p(E) may 
be evaluated exactly from the dynamical rule for general 
restitution coefficient a and coupling strength /. Later, 
we shall consider the generating function of p(E), using a 
procedure which is formally possible for the general case, 
but which will be fully solved only for limiting cases. 

The first moment of p{E) is the average energy, which 
in analogy with granular materials is denoted as the gran- 
ular temperature, Tq = (E). This is evaluated by aver- 
aging Eq. (JJJ over the whole system: 

T G (t + dt) = (l-Tdt)T G (t) + (l-f)TdtaT G (t) 

+ fTdt±(T G (t)+T B ). (2) 



In the steady state T G (t + dt) = T G (t), and therefore 

T G = T B / [2a - 1 + 2(1 - «)//] . (3) 

While the model described above is extremely simple, 
it is interesting to note that it captures, at least qualita- 
tively, some aspects of an actual driven granular gas. In 
Eq. J2J we see that T G is always smaller than Tb and 
depends not only on the dissipation through the resti- 
tution coefficient a, but also on the details of the cou- 
pling to the bath through the coupling strength / (T G 
coincides with Tb only for the two non-dissipative lim- 
iting cases: conservative interactions (a = 1) and inter- 
actions only with the bath (/ = 1)). To see this for 
a real granular gas (in 2D, for convenience), we esti- 
mate the ratio T g /Tb by the following mean- field en- 
ergy balance calculation (see also Q). Consider a gas 
of grains of mass m, diameter D, and restitution coeffi- 
cient e, at volume fraction <E>. The mean time between 
collisions for a grain of energy E is and the 

mean energy dissipated per collision is proportional to 
(1 — e 2 )T G . Consequently, the average energy loss rate 
due to collisions is K<s ^ 1 Zf: ^ Tt! 2 where K is a dimen- 

D^/rn Or 

sionless constant. This is balanced in the steady state by 
gain due to driving by, and loss due to friction with the 
bath, 2 7 (T B - T G ), yielding T B ~T G - AT B 1/2 T G /2 = 0, 
which may be solved for T G /Tb as a function of the sin- 
gle parameter A(e, C) = — e 2 ) (where we use the 

dimensionless coupling strength C = jD^Jlf^). We note 

that results of molecular dynamics simulations agree with 
this estimate, as can be seen in Fig. QJi. This granular gas 
energy balance is more complicated than the calculation 
leading to Eq. J2Jl, since the rate of dissipative interac- 
tions in a granular gas depends not only on the coupling 
strength C and volume fraction $, but is dynamically 
determined by the typical energy T G through the grain 
velocities. Moreover, density and energy correlations in 
granular gases likely break the validity of such a mean 
field approximation, especially at high volume fraction 
and low restitution coefficient, or for higher moments of 
the energy. 

For our model the frequency of dissipative interactions 
is determined by the controllable parameter /, and since 
our model is inherently correlation free, any moment of 
the energy may be exactly calculated by taking the av- 
erage of any power of Eq. (|TJ . This yields the following 
recursion relation for (E n ) in the steady state in terms 
of all lower moments and the known moments of Eb, 
{E n B l ) = mlTg: 

(E n ) = (n + l-f-2a n (l-f))- 1 [f(E n B ) 
+ E ( m ) ( Em ) (f( E B- m ) + (1 - f)a n {E n - m ))}- (4) 

m=l ^ ' 

In order to obtain a full solution of the energy dis- 
tribution p(E), we introduce the generating function 
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FIG. 1: a) The ratio Tq/Tb vs. the scaled parameter 
A(e,&,C), defined in the text. Results from numerical sim- 
ulations with different combinations of the restitution coef- 
ficient e, the volume fraction O and the dimensionless cou- 
pling strength C (symbols) agree with the mean field cal- 
culation (solid line). Results with 400 particles are shown, 
however similar results were obtained with up to 12800 par- 
ticles, b) Normalized energy distributions from simulations 
with $ = 0.1 and C = 0.01 (dashed lines), compared to the 
Boltzmann distribution (thick gray line). 



g(X) = (e~ XE ) = J °° e~ XE p(E)dE. Averaging the expo- 
nential of Eq. Q over the whole system and considering 
the steady state yields 

(e- XE ) = {1-Tdt)(e- XE ) + {1- f)Tdt f (e~ zaXE ) 2 dz 

Jo 

+ fTdt f (e~ zXE )(e- zXEB )dz, (5) 
Jo 



dissipative limiting cases a — 1 and / = 1. Furthermore, 
the solution for the limit of maximal dissipation (a = 0) 
is 5 (A) = 2 Fi(l, 2, 2 - /, -XT B ) ■ (XT B + 1), where 2 F X 
is the Gauss hypergeometric function. The theoretical 
results of these two limiting cases (a — and a — 1) are 
plotted in Fig. |21 together with results for < a < 1 ob- 
tained from a Monte-Carlo numerical simulation of the 
model with 1000 particles. The numerical results indi- 
cate that the theoretical results known for a = and 
a = 1 bound the family of solutions for all intermedi- 
ate values < a < 1. Note that the characteristic form 
of the energy distribution in this exactly solvable model 
qualitatively resembles that of granular gases (see Fig. 
\3p), for which only approximate solutions exist. 




FIG. 2: The generating function (a) and normalized energy 
distribution (b) for coupling strength / = 0.5. Dashed lines 
are results of numerical simulations for various values of the 
restitution coefficient a, and thick gray lines are theoretical 
solutions for the two limiting cases. 



9(A) = (1 - /) f 9 2 (za\)dz + f f g( l A) f , (6) 
Jo Jo zAl B + L 

where we have used the generating function of the ex- 
ponential energy distribution in the bath, <?b(A) = 
e- XEB Tg 1 e- EB / TB dE B = {XT B + l)" 1 . We now 
transform the integral equation JHJ into a differential 
equation by change of variables to zX in both integrals, 
multiplication by A and differentiation by A, yielding 



A 



dg(X) 
d\ 



(1 - f)g 2 (a\) 



f 



AT S + 1 



1 5(A). (7) 



We first note that the generating function of the Boltz- 
mann distribution, g B (X), solves Eq. (JJJ for the two non- 



We will now use FD relations in order to demonstrate 
one observable deviation from equilibrium behaviour in 
dissipative systems due to their non-Boltzmann energy 
distribution. In order to measure FD relations in their 
simplest form we introduce an internal degree of freedom 
per particle, Xi, which is coupled to a uniform exter- 
nal field, F, so that the total energy of every particle is 
Ui = Ei — XiF, where Ei is now referred to as the kinetic 
energy. We consider the relation between the fluctuation 
in x, {Ax 2 ) = (x 2 ) — (x) 2 , and its susceptibility with 
respect to changes in F, %= ^fr-. 

We assume driven dissipative dynamics for the kinetic 
energy together with non-dissipative exchange of kinetic 
energy and internal energy, as described by the following 
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equations of motion, 



Ei(t + dt) 



yah 



Ei(t) 
za{Ei{i) + Ej(t)) 
z(Ei(t)+E B ) 
z{Ei(t) -XjF) 



probability : 
1 - Tdt 
(1 - f)Tdt 
/(l - h)Tdt 
fhTdt 



(8a) 



Xi(t + dt) 



value : 
Xi (t) 



probability : 
1 - fhTdt 
fhTdt 



(8b) 



where h is a parameter introduced to describe the ra- 
tio between interactions with the bath and interactions 
between the two types of degrees of freedom. By calcu- 
lating the first two moments of Eq. (|8b|l one sees that 
X = and (Ax 2 ) = That is, the FD relation, 

(Aa: 2 ) = xTfd, is satisfied with an effective temperature 

IE 2 ) 

Tp £> = j^gy , which probes the distribution of the kinetic 



energy E, determined from Eq. (|8a|l . Note that \ an d 
(Aa: 2 ) diverge as F — > (where Tfd is normally defined), 
however their ratio is finite and independent of F . 

Tfd, which characterizes fluctuations, results from the 
second moment of the energy distribution, and is gener- 
ally different from Tq, which describes the first moment. 
T FD and T G coincide only if (E 2 ) = 2(E) 2 , which is the 
case for the Boltzmann distribution. For dissipative sys- 
tems far from equilibrium Tfd is larger than Tq, and 
their ratio is given for our model in the limit h — > by 



T, 



FD 



T, 



G 



(E 2 ) 
2(E) 2 



[2(2 -/)-2a(4 -/)(!-/) 



+a 2 (4 - 3/)(l - /)]/[/(3 - / - 2a 2 (l /))]. (9) 

In dissipative systems with strong coupling (/ ss 1) and 
large restitution coefficient (awl) the energy distribu- 
tion is close to exponential, hence the values of Tfd and 
Tq are similar (but not identical). Generally, the ratio 
Tfd/Tg may reach any value larger than one, and di- 
verges for vanishing coupling strength (/ — ► 0). As has 
recently been predicted by kinetic theory 0] , we expect 
Tfd to be larger than Tq in granular gases as well, where 
the energy distribution is non-exponential. In the cases 
studied numerically the energy distributions were only 
slightly non-exponential, resulting in small differences be- 
tween Tfd and Tq, which explains their seeming coinci- 
dence. 

In conclusion, we have presented a simple dissipative 
model, solved it in terms of all energy moments in the 
general case, and obtained an exact expression for the 
generating function in the maximally dissipative limit. 
Although our model is inspired by granular gases, we be- 
lieve it may have relevance to a broader class of driven 
dissipative systems. A dissipative system coupled to a 
heat bath is very different from a conservative system 



coupled to the same bath. Not only is the granular tem- 
perature lower than the bath temperature, but it also de- 
pends on the coupling details. In addition, FD relations 
hold with an effective temperature Tfd which character- 
izes the second moment of the energy distribution, and is 
larger than the granular temperature, Tq. These exactly 
coincide in equilibrium, where the energy distribution is 
exponential, however generally differ in dissipative sys- 
tems. We expect that careful numerical studies of granu- 
lar gases will show this heretofore unobserved difference. 
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